???+3 



Price and Prichard: Non-linear Response of the Magnetosphere 



Acknowledgements. We wish to thank J. Theiler, N. Hunter, 

D. Vassiliadis, A. Klimas, S. Sharma and T. Chang for many 
useful discussions. We thank NSSDC for assistance in acquir- 
ing the ISEE-3 data, and acknowledge S. J. Bame (plasma) and 

E. J. Smith (magnetometer) as Principal Investigators on those 
ISEE-3 instruments. One of us (D. P.) would like to acknowl- 
edge the support of a UAF Chancellors Graduate Fellowship 
and an AWU/DOE Thesis Parts Fellowship. This work is 
supported by NSF grant ATM-92I3522. 

References 

Baker, D. N., R. D. Zwickl, S. J. Bame, E. W. Hones, Jr., B. T. 
Tsurutani, E. J. Smith, and S.-I. Akasofu, An ISEE-3 high 
time resolution study of interplanetary parameter correla- 
tions with magnetospheric activity, /. Geophys. Res., 88, 
6230, 1983. 

Baker, D. N., A. J. Klimas, R. L. McPherron and J. Biichner, 
The evolution from weak to strong geomagnetic activity: an 
interetation in terms of deterministic chaos, Geophys. Res. 
Lett, 17,41, 1990. 

Bargatze L. R, D. N. Baker, R. L. McPherron, and E. W. Hones, 
Jr., Magnetospheric impulse response to many levels of ge- 
omagnetic activity, /. Geophys. Res., 90, 6387, 1985. 

Brock, W. A., Nonlinearity and Complex Dynamics in Eco- 
nomics and Finance, in The economy as an evolving complex 
sy.stem, vol V of SFI Studies in the sciences of complexity, 
P. W. Anderson, K. J. Arrow and D. Pines, eds., p. 77, 
Addison- Wesley, 1988. 

Casdagli, M., Nonlinear prediction of chaotic time series, 
Physica D, 35, 335, 1989. 

Casdagli, M., Chaos and Deterministic vs. stochastic nonlin- 
ear modeling, 7. R Stat. Soc. B, 54, 303, 1992a. 

Casdagli, M., A dynamical systems approach to modeling 
input-output systems, in Nonlinear modeling and forecast- 
ing, vol XII of SFI Studies in the sciences of complexity, 
M. Casdagli and S. Eubank, eds. p. 265, Addison- Wesley, 
1992b. 

Clauer, C. R., R. L. McPherron, C. Searls, and M. G. Kivel- 
son. Solar wind control of auroral zone geomagnetic activity, 
Geophys. /?e.s. Lett., 8, 915, 1981. 

Eckmann, J.-P, and D. Ruelle, Fundamental limitations for es- 
timating dimensions and Lyapunov exponents in dynamical 
systems, Physica D, 56, 185, 1992. 

Eubank, S., and J. D. Farmer, An introduction to chaos and ran- 
domness, in ] 989 Lectures in Complex Systems, SFI Studies 
in the Sciences of Complexity, Lect. Vol 2, Erica Jen, ed., p. 
75, Addison- Wesley, 1990. 

Farmer, J. D. and J. J. Sidorowich, Exploiting chaos to predict 
the future and reduce noise in: Evolution, Cognition, and 
Learning, Y. C. Lee, ed., p. 277, World Scientific, Singa- 
pore, 1988. 

Gonzalez, W. D., B. T. Tsurutani, A. L. C. Gonzalez, E. J. 
Smith, F. Tang, and S.-I. Akasofu, Solar wind - magneto- 
sphere coupling during intense magnetic storms, J. Geophys. 
Res., 94, 8835, 1989. 

Hunter, N. P., and J. Theiler, Characterization of nonlinear 
input-output systems using time series analysis in Proceed- 
ings of the 1st Experimental Chaos Conference, S. Vohora, 
M. Spano, M. Schlesinger, L. Pecora and W. Ditto, eds., p. 
54, World Scientific, Singapore, 1992. 

Kamide, Y., and S.-I. Akasofu, Notes on the auroral electroject 
indices. Rev. of Geophys. Space Phys., 21, 1647, 1983. 

KHmas, A. J., D. N. Baker, D. A. Roberts, D. H. Fairfield, 
and J. Biichner, A nonlinear dynamical analogue model of 
geomagnetic activity, 7. Geophys. Res., 97, 12253, 1992. 



Ogilvie, K. W., A. Durney and T. von Rosenvinge, Descrip- 
tions of experimental investigations and instruments for the 
ISEE spacecraft, IEEE Trans. Geosci. Electron., GE-16, 
151, 1978. 

Prichard, D., and C. P. Price, Spurious dimension estimates 
from time series of geomagnetic indices, Geophys. Res. 
Lett., 19, 1623, 1992. 

Prichard, D., An investigation of the non-linear response of 
the magnetosphere, thesis submitted in partial fulfillment of 
requirements for M. S. in Physics at UAF, 1993 (unpub- 
lished). 

Roberts, D. A., D. N. Baker, A. J. Klimas and L. F. Bargatze, 
Indications of low dimensionality in magnetospheric dy- 
namics, Geophys. Res. Lett, 18, 151, 1991. 

Roberts, D. A., Is there a strange attractor in the magneto- 
sphere?, /. Geophys. Res., 96, 16031, 1991. 

Sauer, T., J. A. Yorke and M. Casdagli, Embedology, J. Stat. 
Phys., 65, 579, 1991. 

Shan, L., P. Hansen, C. K. Goertz, and R. A. Smith, Chaotic 
appearance of the AE index, Geophys. Res. Lett., 18, 147, 
1991a. 

Shan, L. H., C. K. Goertz and R. A. Smith, On the embedding- 
dimension analysis of AE and AL time series, Geophys Res. 
Lett., 18, 1647, 1991b. 

Takens, R, "Detecting strange attractors in turbulence", in 
Dynamical systems and turbulence. Vol 898 of Springer lec- 
ture notes in Mathematics, p. 366, Springer- Verlag, Berlin, 
1981. 

Theiler, J., Spurious dimension from correlation algorithms 
applied to limited time series data, Phys. Rev. A., 34, 2427, 
1986. 

Theiler, J., Estimating Fractal dimension, J. Opt. Soc. Am. A, 
7, 1055, 1990. 

Theiler, J., B. Galdrikian, A. Longtin, S. Eubank, and J. D. 
Farmer, Using surrogate data to detect nonlinearity in time 

series, in Nonlinear Modeling and Forecasting, vol XII of 
SFI Studies in the sciences of complexity , ed. M. Casdagli 
and S. Eubank, p. 163, Addison- Wesley, 1992; Theiler, J., 
S. Eubank, A. Longtin, B. Galdrikian, and J. D. Farmer, 
Testing for nonlinearity in time series: the method of surro- 
gate data, Physica D, 58, 77, 1992. 
Tsurutani, B. T., J. A. Slavin, Y. Kamide, R. D. Zwickl, J. H. 
King, and C. T. Russell, Coupling between the solar wind 
and the magnetosphere: CD AW 6, J. Geophys. Res., 90, 
1191, 1985. 

Tsurutani, B. T, M. Sugiura, T. lyemuri, B. E. Goldstein, W. 
D. Gonzalez, S.-I. Akasofu and E. J. Smith, The nonlinear 
response of AE to the IMF driver: a spectral break at 5 
hours, Geophys. Res. Lett., 17, 279, 1990. 

Vassiliadis, D. V., A. S. Sharma, T. E. Eastman, and K. Pa- 
padopoulos. Low-dimensional chaos in magnetospheric ac- 
tivity from AE time series, Geophys. Res. Lett., 17, 1841, 
1990. 

Vassiliadis, D. V., A. S. Sharma, T. E. Eastman, and K. Pa- 
padopoulos, Lyapunov exponent of magnetospheric activity 
from AL time series, Geophys. Res. Lett, 18, 1643, 1991. 

C. P. Price and D. Prichard, Physics Department, University 
of Alaska Fairbanks, AK 99775 

(Received February 1, 1993; 
accepted March 23, 1993.) 



Price and Prichard: Non-linear Response of the Magnetosphere 



??? 



N=1440 min., x=15 min., W=50 min., D|,=10 




r(nT) 

Fig. 3. The BDS statistic as a function of radius for the AE 
data set (solid line) and for the surrogates (dashed lines). The 
delay time is t = 15 minutes, the value of W is 50 minutes, 
and the embedding dimension rn = 10. 

Figure 4 shows the singular value decomposition dimension 
estimate i?^^^ as a function of the number of neighbors for the 
AE data set (solid line) and for the surrogates (points). The 
delay time r is 1 minute, the embedding dimension is = 
15 and the least significant eigenvalue is chosen arbitrarily to 
be 50 times smaller than the largest eigenvalue. The results 
are not affected significantly by variations in the embedding 
dimension (we have calculated out to = 60) or the least 
eigenvalue size. We caution the reader from iriferring from 
this figure that the AE data set shows a dimension of roughly 7 
(or that the stochastic surrogates show a dimension of roughly 
10!). Again, we are using the singular value dimension esti- 
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Fig. 4. The singular value decomposition dimension estimate 
versus the number of neighbors for the AE data set (solid 
line) and for the surrogates (points). The delay time is 1 
minute, the embedding dimension is = 15, the value of W 
is 50 minutes, and the least significant eigenvalue is 50 times 
smaller than the largest eigenvalue. 
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Fig. 5. Normalized one step ahead prediction error (Ej) 
versus the number of neighbors. The delay time is r = 15 
minutes, and the embedding dimension is = 8 for both 
the input and output data sets. The training interval for the 
nonlinear predictor contains 30000 data points. 

mate as a non-linear statistic, not as a measure of the effective 
degrees of freedom. What Figure 4 does show is a statistically 
significant difference between the AE data set and the ensemble 
of surrogates. 

In Figure 5 we show the normalized prediction error versus 
the number of neighbors used to compute the local maps. The 
delay time r is 15 minutes and the embedding dimensions 
are D^^ = D^2 — 8- The normalized error shows a distinct 
minimum of = 0.137 at roughly 600 neighbors (out of a 
30000 point training set). For very large numbers of neighbors 
(i.e. a linear predictor) the normalized error approaches E^ = 
0.15. Therefore, the non-linear predictions are roughly 10% 
better than linear predictions. 

Conclusions 

On the day 30 October 1978, the IMF had a nearly constant 
southward value oi Bz = - 10 nT, and the magnetospheric ac- 
tivity (as indicated by the AE index) was quite high. Although 
the magnetosphere is not an autonomous system, when the in- 
put forcing function is relatively steady, the system may have 
time to converge to an attractor. Hence, while in general the 
magnetosphere must be analysed as a non-linear input/output 
system, for this day (when the input function is approximately 
a point in state space) we can use both the standard dimensional 
estimates and non-linear prediction methods. 

Our approach has been solely to determine whether or not 
the magnetosphere gave evidence of a deterministic non-linear 
response on that day. We are not yet attempting to quantify the 
non-linearity. Using convenient non-linear statistics, we have 
compared the AE data set to an ensemble of surrogate data 
sets. In the same vein, the non-linear prediction analysis is a 
comparison of the non-linear prediction to a linear prediction 
of the output AE signal from the input signal vBg. 

We have found diat the non-linear statistics show signifi- 
cant differences between the AE data set and the ensemble of 
surrogates. Further, we find that the non-linear predictor is 
somewhat better than a linear predictor. The combination of 
these results leads us to conclude that, for 30 October 1978, 
there is evidence for deterministic non-linear coupling between 
the solar wind and the magnetosphere. Further work is needed 
to identify the nature and characteristics of that coupling. 
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Fig. 1. The time series for the solar wind input vBs and 
magnetospheric response AE for 30 October 1978. 

the autocorrelation time. C should scale (over a limited range) 
as: C(r) ~ r'^ where d is the dimension. Takens [1981] sug- 
gested a direct dimensional estimate based on the M distances 
less than some upper cutoff tq. 
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We choose tq to be roughly | the rms deviation of the time 
series; the Takens estimate can be calculated directly from 
the correlation integral [Theiler, 1990]. Another non-linear 
statistic closely related to the correlation integral is the BDS 
statistic: BDS = VN[C{r; m) - C{r\ 1)™] where C{r\ m) is 
the correlation integral evaluated at radius r in m embedding 
dimensions [Brock, 1988]. For stochastic time series, C(r; m) 
= [C(r; 1)]™, so the deviation of \he BDS statistic from zero 
gives a measure of the determinacy of the time series. 

Another dimension estimator is local singular value decom- 
position [Broomhead et al, 1987; Passamante et al, 1989]. 
The nearest n neighbors to a reference point form a hyperel- 
lipsoid with m significant principal axes; iii is the local dimen- 
sion of the attractor manifold. Averaging m over the reference 
points gives the dimension estimate D , ,. j. The choice of delay 
is almost arbitrary. The nx dp^ trajectory matrix 
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\Sxn 



8x2+r 
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is used to form the covariance matrix Z — X'X. Here 8xj is 
the difference between the point Xj and the reference point, m 
is estimated by the number of significant eigenvalues of Z. 

The best approach is to consider the magnetosphere as a 
non-linear input/output system, that is, to simultaneously ex- 
amine both solar wind parameters and geomagnetic indices, 
and compare the output to a non-linear predictor [Farmer and 
Sidorowich, 1988; Casdagli, 1989; Eubank and Farmer, 1990; 
Casdagli, 1992a,b; Hunter and Theiler, 1992]. For a driven 
system with input u{t) and output y{t) we make embedding 



vectors c{t) = [u{t). u{t - r) . . . v{t -{Dp^- 1)-). y{t), y{t- 
-)... ?/(/ - {Dfr,2 - We use the first 29 days of October 

1978 to "train" the phase space \e}. To predict y{l. + 1) we 
first find the k nearest neighbors to e{t) in the state space, then 
we determine where these points go after 1 time step, and then 
find a least squares solution for y{t + 1). 

To determine the quality of the prediction, we randomly 
choose 50 points in the state space during 30 October 1978, 
and predict the value 1 step ahead. After all 50 predictions are 
made we compute the normalized error Ej^i 



1 ^ir^ predicted Vreal'x'^ 



with a the standard deviation of the time series. As k ap- 
proaches the "training" set length this method becomes a linear 
predictor. If there is a non-linear relationship between input 
and output then there will be a minimum of £y for k small 
compared to the length of the "training" set [Casdagli, 1992a]. 
Results 

We have examined non-linear statistics for the 30 October 
1978 AE data set and for ten surrogates to that set. Figure 
2 shows the dependence on the Takens dimension estimate 
Dxaiens ^ functlou of thc embedding dimension D^ for 
the AE data set (solid line) and for the surrogates (points). 
The delay time r is 15 minutes, and W is 50 minutes. We 
caution the reader from drawing the inference from this figure 
that the AE data set shows a dimension of roughly 3; in fact, 
the slopes of the correlation integrals are neither convergent as 
Dp is varied nor flat at any radius. We remind the reader that 
we are using the Takens dimension estimate as a non-linear 
statistic, not as a measure of the effective degrees of freedom. 
What Figure 2 does show is a statistically significant difference 
between the AE data set and the ensemble of surrogates. 

Figure 3 shows the BDS statistic as a function of radius 
for the AE data set (solid line) and for the surrogates (dashed 
lines). The delay time r is 15 minutes, W is 50 minutes, 
and the embedding dimension m = 10. Figure 3 also shows a 
statistically significant difference between the AE data set and 
the ensemble of surrogates. 
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J, versus the 

embedding dimension Dp for the AE data set (solid line) and 
for the surrogates (points). The delay time is r = 15 minutes, 
and the value of W is 50 minutes. 
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Abstract. Previous efforts to find evidence of deterministic 
nonlinear dynamics in the global geomagnetic system have 
treated the geomagnetic system as autonomous. However, the 
geomagnetic system is strongly driven by the stochastic so- 
lar wind. We consider the response of the magnetosphere, as 
given by the AE index, for one day when the IMF had a nearly 
constant southward value. Using both a series of non-linear 
statistics and non-Unear prediction of the response to the input 
signal vBs,we find that there is some evidence for determin- 
istic non-linear response of the Earth's magnetosphere on that 
day. 

Introduction 

There have been numerous studies [e.g. Clauer et al, 1981; 
Baker et al, 1983; Bargatze et al, 1985; Tsurutani et al, 
1985, 1990; Gonzalez et al, 1989; Vassiliadis et al, 1992b] 
trying to relate solar wind parameters to geophysical indices 
(e.g. AE, AL, Kp, etc.) These studies show that the mag- 
netospheric system is neither periodic nor quasiperiodic, and 
that the magnetosphere acts as a type of low-pass filter on the 
solar wind input. This suggests the possibility of a non-linear 
magnetospheric response. These studies also show that only 
part of the AE variation can be explained as being a direct 
linear response to the solar wind. 

Recently there has been much interest in determining 
whether the Earth's magnetospheric system behaves like a low 
dimensional (i.e. deterministic) non-linear system [Vassiliadis 
etal, 1990; Shwetal, 1991a,b; Roberts era/., 1991; Roberts, 
1991; Vassiliadis efaZ., 1991; Prichard and Price, 1992]. These 
studies, which focus on calculating the properties of a possible 
strange attractor in the Earth's magnetosphere, do not yield 
clear evidence for deterministic non-linear behavior, although 
some non-linear circuit analogue models show promising re- 
sults [Baker et al, 1990; Klimas et al, 1992]. 

All of the previous non-linear studies of the magnetospheric 
response have focused solely on the the state of the magneto- 
spheric system. There are problems with this approach. First 
and foremost, calculations intended to discover the properties 
of a strange attractor (viz. the fractal dimension) assume tac- 
itly that the system is either autonomous or periodically driven. 
When a system (such as the Earth's magnetosphere) is driven 
stochastically, the trajectory in phase space is usually not able 
to converge to an attractor. Thus, attempts to determine the 
non-linear properties of the system solely from the response 
of the system (no matter how complete the diagnostics are) 
will not be meaningful, since the strange attractor is in general 
never delineated. Clearly, it is imperative to consider simul- 
taneously the solar wind driver, and the only way to do so in 
general is to use non-linear input/output prediction techniques 
[Farmer and Sidorowich, 1988; Casdagli, 1989; Eubank and 
Farmer, 1990; Casdagli, 1992a,b; Hunter and Theiler, 1992]. 
Second, since the results of non-linear calculations of any kind 
on real data sets frequently give less than clear-cut results, it 
behooves one to have a standard of comparison, for which we 
will use surrogate data sets [Theiler et al, 1992]. 

In this paper, we analyse the response of the magnetosphere, 
as given by the AE index, for the day 30 October 1978, when 
the interplanetary magnetic field (IMF) had a nearly constant 
southward value. Our choice is motivated by the following 
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reason: only for periods when the forcing function is relatively 
steady, i.e. the input function reduces to a point in state space, 
is it possible for the trajectories in phase space to have time to 
converge onto a strange attractor. Thus, this interval may be 
seen as a bridge between the previous analyses and non-linear 
prediction analyses. We present the results of two distinct 
analyses of the magnetospheric response: first, we make the 
standard dimensional estimates (via correlation integral and 
local singular value decomposition), and second, we make a 
non-linear prediction of the response to the input signal vB.,. 
Data Sets and Methods of Analysis 

The response of the magnetosphere is taken to be given 
by the auroral electrojet index AE, with one minute resolution. 
Although there are problems with this choice [see e.g. Kamide 
and Akasofu, 1983; Prichard, 1993], the use of the electrojet 
indices has been standard for this type of investigation. Here, 
the solar wind input is taken to be the product of the southward 
component of the IMF (Eg — \[\Bz\- Bz\) and the solar wind 
speed, calculated from ISEE-3 measurements [Ogilvie et al, 
1978] with one minute resolution. Figure 1 shows the solsir 
wind input vB., and the magnetospheric response AE time se- 
ries for 30 October 1978; note during a large part of the day 
the IMF had a nearly constant value of = - 10 nT. For the 
non-linear prediction, we will use the data of the entire month 
October 1978. Although the 30 October 1978 AE data set 
length is rather short (1440 elements) for dimension estima- 
tion [Eckmann and Ruelle, 1992] we are using the dimension 
estimates as convenient non-linear statistics (see below). 

We have also constructed surrogate data sets (nondetermin- 
istic sets with the same statistical properties, viz. power spec- 
trum, autocorrelation, etc.) for the AE data in order to provide 
a standard for comparison. Surrogates can easily be made 
by taking the Fourier transform, randomizing the phases, and 
then inverting the transform. Our surrogates also have the 
same amplitude distribution in the time domain as the original 
data [Theiler et al. 1992]. We employ the surrogate data sets 
to determine the significance of the results of our analyses. If 
the results of the non-linear statistics (here, dimension esti- 
mates) for the AE data set vary meaningfully from the results 
for the surrogate data sets, one may assume the presence of a 
deterministic non-linear process in the magnetosphere. 

We assume that the signal {e.g. AE) is due to a non-linear 
dynamical process involving an unknown number of variables. 
Using the method of time delays (Packard etal, 1980; Takens, 
1981) we reconstruct the attractor by creating a set of 
dimensional vectors from the original time series {x[i)Y. x^ — 
[x{i), x{i + T),x{i + 2T), . . . , x{i + {D ^ - l)r)] with t the time 
delay, and the embedding dimension. This reconstruction 
will have the same topological and dynamical properties as the 
original attractor as long as > 2D, where D is the attractor 
dimension (Saueref a/., 1991.) 

One estimate of the attractor dimension uses the scaling of 
the number of points inside a ball as the radius of the ball (r) 
increases. One computes the correlation integral [Grassberger 
and Procaccia, 1983]: 

- [n-w)'n-w^,)% g - ^'^-^ - ^^^11) 

Removing W points on either side of the reference point elim- 
inates autocorrelation effects [Theiler, 1986]; for = 1 we 
recover the standard correlation integral. W is on the order of 
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